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Abstract 


Wc report on the results of the first year's activities under NASA grant 
5254, covering a period from March 15th 1978 to March 14th 1979, "Stochastic 
modeling of the time-averaged equations for climate dynamics". We discuss two 
analyses, both based on a simplified set of barotroplc equations for representing 
large scale non-linear atmospheric circulation characteristics. 

We have essentially completed study of the statistical properties of this 
set of equations for small values of a parameter , corresponding to the 
physical case of large zonal westerlies and relatively weak eddy motions. We 
have also started an investigation of the effect of seasonal type forcing on the 
solution of these equations, and some preliminary numerical results are presented. 

In addition to describing the formal analyses, we Include a more general 
discussion on the significance of the results In the context of the use of the 
simplified equations, and indicate the directions along which we expect progress 


to be made. 


Introduction 


The first year's work performed under NASA grant 5254 has concentrated on 
study of the statistical properties of a particularly simple model for large scale 
atmospheric motions - a three wave component barotropic system, first introduced 
in the meteorological literature by toren?'^^, the so called "minimum hydrodynamic 
system" of equations. The immediate objective of our program was to investigate 
the long term statistical prop<^rtles of this model, employing time-averaged 
quantities as the dependent state variables. Introduction of the latter, it was 
hoped, could lead to direct means for obtaining detailed climatic states more 
efficiently than by computing average statistics from long term integration of 
deterministic general circulation models. 

Although such a basic approach to climate is appealing on fundamental grounds, 
it leads to considerable analytic complexity because of the multiplication in 
the number"of differential equations that result from the time-averaging procedure; 
we are in fact confronted with the need to deal with a formally infinite number of 
coupled non-linear differential equations defining the statistical moments of the 
system. It was expected that the problem of closure of this set could be tackled 
more easily by using the simplified minimum equations as a model for the more 
complete and realistic system, not just because of ease of computation of numerical 
solutions, but primarily because analytic solutions for the minimum equations are 
available, and these can be used as a check on approximate methods of closure of 
the moment equation set. 

One approximation method that held promise (based in part on experience gained 
with the same system of equations applied to weather forecasting problems' ') was 
a linearisation of the moment equations and truncation at a fairly low order for 
estimation of small climatic changes. A large part of the first year's work 
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stemmed from this approach, but, as we shall describe shortly, during the course 
of the Investigations an Important property of the minimum equations was dis- 
covered that largely eliminates the need to consider a small climatic change 
linearisation approximation for closure. 

A second theme. Introduction of seasonal variation Into the study of climate 
statistics, was also started. Here again we have used the minimum equation system 
as a model and simulated seasonality by Introduction of additional forcing terms 
In the equations. 

CloS’are of the moment equations for the minimum equation s 

In the process of calculating magnitudes of higher order moments of the 
minimum equation solutions, an unexpected characteristic was discovered, namely 
a systematic reduction In magnitude of the moments with Increasing order, thus 
enabling us to formulate a closure technique for the Infinite set of moment equa- 
tions that we can demonstrate gives good approximations to exact solutions. This 
result applies for small values of the parameter (to be defined below) 

corresponding to the case of relatively large Westerlies and low eddy energy, 
matching conditions that Loren?, originally suggested best approximated mid- 
latitude conditions on the earth. We have completed a comprehensive analysis of 
the reasons for this useful result, and present next a summary of the approach and 
the analytic results. 

The equations governing the temporal development of the amplitudes of the 
three largest scale Fourier components of the barotropic flow in a rectangular 
domain have been shown by Lorenz to be 
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where 
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Cj - (1 - o^)/2m 

and a - k/l » the ratio of the meridional to zonal dimensions of the rectangular 
domain. We note that the important Coriolis grandient effects are omitted in 
this model. The time averaged form of the moment equations can be shown to be 
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where the moments are defined as 
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As shown by Lorenz^^^ exact analytic solutions for (1) are given by 

ft* » ft? Sh /A/f.f*), »<.) •<>' (5) 
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where dn , sn , cn are elliptic functions, modulus , and h , t* and 
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are given in terms of the maximum amplitudes ft, , ft^, ftj , by 
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The parameter k^ determines the periodicity of the elliptic function 

solutions, which in turn corresponds to a vacillation phenomenon, analogous 

to variations in zonal index that occur in the real atmosphere. Larger k^ 

values corresponds to slower cycling, the periodicity becoming Infinite as k 

(1) 2 

approaches 1. Lorenz chose a value of k^ « 0.2 as representative of typical 
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mid-latitude conditions; this yields a vacillation period of 5.6 days and, with 
the parameter "t- * 2 , a sonal wave length for the Fourier component of 5000 km. 

Exact analytic results can be calculated from (1) for moments, defined over 
an averaging time T . We list values up to the third order moments: 
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where K and £ are complete elliptic Integrals of the the first and second kinds. 

If we let the averaging time approach infinity, a number of the moments defined 
by (7) vanish, and the remainder are given by 
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In this case the moment equations (2) reduce to an algebraic set 
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and we should be able to derive the solutions (8) from (9) , once the two 

Integral invariants of^the system - energy and enstrophy - are prescribed. 

The amplitudes A* . A* and A* are in fact calculable from the energy V 
and enstrophy E : 
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We have found no general technique for solving (9)» given E and V, and 
indeed the situation is made more perplexing by the fact that some of the 
original differential equations (2) are no longer Independent in the steady state 
form, (9). For example, the three time-dependent equations for the second order 
moments 
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reduce to the single equation 
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when the time derivatives vanish. Thus far we see no profound reason for 

studying this anachronism and have not tried to understand structure of the 

steady state system. As we appreciate the subtlles of the statistical behavior 

of the minimum equations we may wish to return to the problem. 
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For small values of , we can expand the elliptic function solutions of 
equation (8) (the Infinite averaging time case) , yielding 
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and from studying tha order of magnitudes (In k^) of the terms In (13) we see 
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that the first and second order moments are In fact given to order by 
solution of the closed set of equations 
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Our numerical work with the minimum equations has centered on choice of 

kp ** 0.2 , as was the origlml study by Lorenz. The approximation for small 

k that we are discussing Is quite adequate to cover this situation, as can 
o 

be seen In the graphs plotted In figure (1). Here we show approximate solutions 
for some of the moments corresponding to equation (14) . We give the variation 
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with of the firet order approximation for and and the 

second order approximation for , 

If we go on to look at the structure of the moment solution system for 
finite T values, again with small values of , we find a somewhat more 
complex ordering, but still of a nature that enables us tc close the moment 
equation set. For example, from the analytic forms (7), using the expansions 
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A systematic procedure for developing the approximation to the moment 

2 

equations for small Kq can be devised by considering an expansion of the 

2 

dependent variables as power series in Rq . The complete development is rather 
laborious, and we shall only outline a few results. We first simplify the algebra 
by introducing non-dimensional variables defined as 
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Wien power series expansions for the moments are Introduced into (21) , a 
series of closed coupled equations result. He cite only the lowest order result. 


12 


Thus, with the expansions 
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which is solvable as 
We note that to 


a closed set, given imtinl values lor the moments 
this order, we have no values, lor example for the 
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covariances i well as for several of the higher order moments. 

To obtain estimates of the former, we need to go to the next higher order equations 


namely 
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in which the F* represent complicated functions of the indicated variables which 
we have no need to specify here. Equation (24) forms another closed set for a 
number of the higher order mometits, as well as the two covariances oj^and ffj’jthat 
were not obtained by solution of equation (23) . 

Earlier we presented some results on the error that results from the small 
assumption (figure (1)) for the infinite time averaged case. Finite time 
averaging does not worsen the situation, as we illustrate in figure (2) , which 
plots that variation in time (in units of the periodicity IK of the elliptic 
function) of the error in ^,^‘*^for various values of the averaging time, T . 


14 




1 


The fractional error rarely exceeds 2%. 

The above Illustrates the systematic procedure that can be used to develop 

solutions for the moment equations, ^^lether the expansion (22) is appropriate 

for equations more general than the minimum equations is not yet knovm, though 

2 

its basis in the small k* assumption appears to have a logical physical foundation 
If Indeed extendable, the result is of great importance, providing at first sight 
a means for dealing with the very formidable problem of closing the moment equation 
system that defines climate means and statistics. 

Seasonal cycle simulation 

The minimum equation system (1) is based on the use of the two-dimensional 
barotropic vorticity equations without forcing or dissipation: 
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where is the stream function. He Introduce a forcing term ■ 

(27) 
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where JZ corresponds to a seasonal cycle, on the right hand side of equation (25). 

Fourier transformation of (25), when so modified, and retention of only 
three Fourier components, as with the development of the unforced minimum equations, 
yields the new set of equations: 
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for the Fourier components of the transformed vortlclty, , With F(t) « 0 

the normal form of the minimum equations (1) Is regained. 

In contrast to the case F(t) ■ 0 , analytic solution of (28) Is not possible 
and numerical techniques have to be used. A numerical Taylor series expansion 
technique for solving (28) has been carried out for a range of values of the 
amplitude of forcing . Stability of the systems seems to be quite sensitive 

to the choice of , but has not been investigated in detail. However, a choice of 

=*6 “2 

F^» 4.8x10 sec gives good stability as well as a seasonal amplitude variation 
that looks reasonable. Figure 4 is an example of some of the results. We show 
the variation during a year of the energy in the system and the behavior of the 
component A 2 (equation 28) . The band of values for A 2 correspond to the limits 
of its oscillation amplitude at a high frequency time scale, corresponding to the 
unforced elliptic function behavior (equations). This frequency now varies during 
the season, In accordance with the formula (30) below; a total of approximately 
ISO cycles occur in one year. The frequency variation is in fact quite slow, 
according to the numerical results, as is the variation in amplitudes of the 
three components Aj^ , A 2 and A^ . To within the accuracy of the calculation, 
the system appears to return exactly to its initial state after a one year cycle. 

At least for the magnitude of forcing that we have used here, the effect of slow 
forcing appears to be rather simple and can be represented in quasi-statlonary 
fashion by introducing appropriate values for the constants that appear in the 
unforced elliptic function solution (5) which we let to vary with the season. 
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From this we draw the Important conclusion that imposition of forcing in our case 
docs not lead to a major change In the spectral distribution of the solution, apart 
from the superposition of the slow amplitude variation for A 2 • (i'here Is no 
corresponding slow change In or A^). In particular, there Is no additional 
long period non-linear oscillation Introduced that could Interact with the forcing 
frequency. 

An empirically determined form for the variation In amplitude and frequency 
of the elliptic function solutions Induced by the forcing we find to be 
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in which the modulus of the elliptic function Is given by 
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where Is the value with no forcing (F(t) - 0) , and 
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It seems that these results should be derivable' . analytically as 
approximations for small values of . We have yet to show that 

this Is the case. . 

Linearisation of the moment equations 

Linearisation of averaged atmospheric state variables for the study of 
climate change should be a valid approximation for small perturbations of climate. 
Note that linearisation for a deterministic (non-averaged) computation of atmos- 
pheric dynamics Is not valid for estimating long term (climatic) changes, because 
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of the large day-to-day variability of the state variables. 

( 2 ) 

Our previous work on stochastic dynamic weather prediction on a linearllied 
approach to assessing the growth in error of weather forecasting models, suggested 
that an analagously useful result might hold in the consideration of small climatic 
changes. Thus we hoped that the statistical moments for small perturbations in 
climate would progressively decrease In magnitude as their order increased, thus 
leading to the possibility of a systematic closure procedure for the corresponding 
moment equations. 

We have been carrying through a check on the validity of this idea using the 
minimum equations as a test case. We have developed numerical algorithms for 
truncation of linearised forms of the moment equations for the minimum equation 
set at the second and third orders. The accuracy of such approximations can be 
judged by comparison with exact solutions for the moments obtainable from the known 
analytic solutions of the minimum equations (equation 5) . 

Persistant programming errors have thus far prevented us for assessing 
whether this approach to closure of the moment equations Is or is not possible. 

Future work 

Although some questions remain concerning the significance of the small k^ 
approximation discussed above, we believe we have now covered this particular 
property of the minimum equations adequately. We have yet to ascertain whether 
the important simplification in the stochastic formulation that it implies is 
extendible to other more complete and realistic climate simulation models. 

A serious concern in planning our future studies relates to the relevance 
to the real world situation of the minimum equation behavior in their original or 

. i 

in modified form (modified by the introduction of a p effect, deterministic or 
stochastic forcing and dissipation). It is clear that the minimum equations 
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are vastly over simplified In many respects, and their utility os partial simulators 
of real atmospheric effects depends a great deal of the mathematical richness of 
their non-llnear properties. Let us briefly discuss two major apparent weaknesses. 
Firstly, we note that all three solutions (5) are phase coherent; but Introduction 

( 3 \ 

of the Coriolis gradient effect, as demonstrated by Pnegle and Robl eliminates 
this oversimple property without destroying analytic simplicity. Secondly, the 
solutions of the minimum equations are strictly periodic, as therefore are also 
their climatic means. This Is quite different from more complicated simulation 

(A) 

models, such as that of Lorenz' as well as the real atmosphere. Non-perlodlclty 
is a fundamental characteristic for long term behavior of non-llnear systems that 
it seems Important for us to duplicate in any simplified models. 

We had, until recently, considered the latter deficiency of the minimum equa- 
tions to be serious enough to cause us to drop their use as simulators of even the 

(5) 

most simple of climatic phenomena. However, a new paper by Dutton has demon- 
strated the radical effect on behavior of the minimum equations when stochastic 
forcing is added. In particular, he shows that the revised system has bifurcation 
properties and multiple solutions, with a degree of complexity quite surprising 
in view of the simplicity of the behavior of the original, unforced system. In 
view of this result we have revised our opinion, and we now propose to continue 
study of the minimum equations. Introducing stochastic forcing as we had originally 
planned, and now also studied by Dutton. In contrast to his work, however, which 
was a deterministic analysis of the time evolution of a single trajectory in phase 
space, we shall look at the stochastic properties of the forced system. This work 
has yet to be started, but wc hope to do so in the coming year. 

With the completion of the studies on seasonal forcing discussed above, and 
the proposed introduction of stochastic forcing, there only remains the introduction 
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of the effect to round out a complete picture of what the three-component 
minimum equatloiia system can yield. It la our present Intent to try to carry 
through such a program to fruition. 

We should also mention here yet another property of the minimum equations 
that has In the past led to criticism of their use as a real world simulator, but 

which recent developments in theory of large scale atmospheric simulation down- 
plays. The minimum equations describe an inherently two dimensional barotroplc 
model for long term atmospheric behavior, thus omitting all the baroclinlc effects 
that are essential in an understanding of the overall atmospimric circulation 
pattern and its energetics, We now see appearing, however, a growing body of 

theoretical and numerical work^^^ suggesting that large scale atmospheric 

(os 

structures, not only on earth but even on Jupiter and Venus can very well be 
simulated in barotroplc lodels. Forcing of such models must, of course, arise 
from baroclinlc effects, as indeed would be represented in our proposed further 
work with the minimum equations, by stochastic forcing terms. 


20 


References 


(1) E, N. Lorenz, Maximum Simplification of the Dynamic Equations. Tollus 12 
243 1960. 

(2) J. A. Laurmann. A Small Perturbation Approximation for Stochastic Dynamic 
Weather Prediction, Xellus 30 404 1978. 

(3) J. Paogle and E. Robl. Ihe Time Behavior of the Probability Density Function 
of Some Simplified Atmospheric Flows. J. Atm. Scl. 34 979 1977. 

(4) E. N. Lorenz. Deterministic Nonperiodic Flow. J. Atm. Sci. 20 130 1962. 

(5) J. G, Vlckroy and J. A. Dutton. Bifurcation and Catastrophe In a Simple, 
Forced, Dissipative Quasi-Geostrophlc Flow. J. Atm. Scl. 36 42 79. 

(6) P. B, Rhlnes. Waves and Turbulence on a Beta-plane. J. Fluid. Mcch. 69 
417 1975. 

(7) G. P. Williams. Planetary Circulations: 1. Barotroplc Representation of 

Jovian and Terrestrial Turbulence, J. Atm, Scl. 35 1399 1978. 

(8) R. Salmon. IVo-Laycr Quasi-Geostrephic TurbwXrnce in a Simple Special Case. 
Geophys. Astrophys. Fluid Dyn. 10 25 1978. 

(9) W. B. Rossow and G. P. Williams. Large Scale Motion in the Venus Stratospher 
J. Atm. Sci. 36 377 79. 


